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(§4) Proc6d6 de reconstruction multi-6chelle de I'image de la structure d'un corps. 



(57) Pour reconstruire la structure d'un corps au moyen 
oSlgorithmes iteratifs on choisit, pour les premieres itera- 
tions de les effectuer avec une resolution moins grande de 
manidre a limiter le nombre de calculs. On montre alors 
que, pour une reduction par trois, voire par quatre de la du- 
r6e du calcul on salt reconstruire I'image. On montre m§me 
que Timage ainsi reconstruite est de meilleure quality. 
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PROCEDE DE RECONSTRUCTION MULTI-ECHELLE DE 
L 0 IMAGE DE LA STRUCTURE D°UN CORPS 

La presente invention a pour objet un procede de 
5 reconstruction multi-echelle de l" image de la structure 
d Q un corps o Elle trouve plus particulierement son 
application dans le domaine medical ou on procede a la 
reconstruction des structures internes de patients sous 
examen, Elle peut trouver neanmoins son application 

10 dans d'autres domaines, notamment dans le controle non 
destructif industriel dans lequel des examens du meme 
type que les examens medlcaux sont pratiques «, La methode 
de 1° invention vise a accroitre la rapidity de la 
reconstruction a trois dimensions quand les algorithmes 

15 de reconstruction sont des algorithmes par iteration „ 

On connait des protocoles d * acquisition utilises en 
medecine pour connaitre les structures internes d'un 
corps o Ces protocoles sont essentiellement les 
protocoles de types tomodensitometrie, RMN, ultrasons, 

20 scintigraphies, Dans ces protocoles on soumet des parties 
du corps d'un patient a un examen*, De maniere par 
exemple a pouvoir reconstruire 1 9 image en trois 
dimensions de la structure interne du corps, un tel 
examen comporte des series de mesures. Chaque serie est 

2 5 distincte d'une autre par la modification d'un parametre 
de la serie de mesures ou experimentation - Par exemple 
en tomodensitometrie, le parametre concerne 
1 8 orientation du tomodensitometreo En RMN le parametre 
peut concerner la valeur d'une impulsion de gradient de 

30 selection de coupe « Dans tous les cas on acquiert ainsi 
plusieurs series de resultats de mesures. Chaque serie 
est representative du corps et de la configuration 
correspondante de la caracteristique d 1 experimentation « 
Pour reconstruire, a partir des series de mesure, 
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1° image a trois dimensions de structures internes sous 
examen, il existe essentiellement deux types 
d'algorithme de reconstruction en imagerie par rayons X« 
Un premier type concerne un calcul de r£troprojection et 
5 filtrage ou encore de reconstruction par trans formee de 
Fourier a plusieurs dimensions o Un deuxieme type, celui 
qui est concerne par 1° invention concerne les methodes 
iteratives de reconstruction dites encore algebriques* 
Ce deuxieme type trouve une application particuliere en 

10 tomodensitom6trie parce que la mise en oeuvre des 
algorithmes du premier type conduit a des solutions 
instables quand il s'agit de reconstruction en volume a 
partir d B un faible nombre de vueso lie principe d'un tel 
algorithme algebrique, dont un perfect ionnement a ete 

15 decrit dans une demande de brevet fran^ais du demandeur 
deposde le 20 janvier 1989 sous le numdro 89 0O676, 
est le suivantp 

On commence d 1 abord par estimer , a priori , 
1° allure de la structure du corps sous examen » Ceci 

20 signifie que, pour tous les elements de volume, ou 
voxels, de ce corps sous examen, on considfere que la 
valeur de la propriety physique testee par 
1 0 experimentation a une valeur connue, par exemple 0, 
Puis on soumet cet objet theorique (ne ressemblant bien 

25 entendu pas du tout a 1 ? objet que l B on cherche a 
reconnaitre) mathematiquement a une serie 

d 5 experimentations semblables a celles auxquelles a 
reellement ete soumis le corps a examiner* La simulation 
ainsi executee conduit elle aussi a des series de 

3 0 "mesures", appelees par exemple repro j ections , de 
1° objet theorique (faux) selon de memes configurations 
de la caracteristique d 8 experimentation <> On compare 
alors serie de mesures a serie de reprojections et, dans 
les series respectives, chaque valeur deux a deux entre 
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elles, de manifere a deduire, a partir de I'ecart 
constate, une modification correspondante de la valeur 
de la grandeur physique dans le voxel concerned On 
modifie en consequence 1* estimation theorique de cet 
5 objeto 

Ceci 6tant fait, on reitere cette operation de 
simulation d g experimentation, et on reitere egalement 
1" operation de comparaison des series de reprojection 
nouvellement obtenues avec les memes series de mesures 

10 r^ellement mesurees dans le corps „ On en deduit un 
deuxifeme ecart avec lequel on remodifie la connaissance 
qu'on avait prealablement obtenu de la structure du 
corps o Ainsi de suite on affine 1* estimation de 1 1 image 
a trois dimensions de la structure du corps o 

15 on doit indiquer que les mesures reelles sont 

acquises avec une resolution donnee et que la 
reconstruction de 1' image du corps sera effectuee selon 
une resolution donnee Egalement o De maniere a bien 
differencier la resolution des mesures et la resolution 

20 de la reconstruction on parlera pour cette derniere, 
dans la suite de cet expose, d°une finesse de 
reconstruction « On sait que la resolution des mesures, 
au moment des acquisitions conduit en correspondance a 
une finesse significative de la reconstruction du corps 

25 sous exameno Une finesse est significative si elle 
correspond ef fectivement a la connaissance la plus 
precise possible, compte tenu de la resolution des 
mesures effectuees. En effet, si on connait 1 8 image 
d'une structure selon une finesse donnee et que, par des 

30 procedes d'agrandissement quelconques, on exprime 
l ff image de cette structure selon une finesse plus 
grande, on n'augmente pas sa finesse significative o En 
effet on ne fait, en definitive, qu 9 utiliser les 
informations acquises. Ceci est assimilable au grain 
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d 11 une photographie en positif qui ne peut jamais avoir 
une resolution superieure au grain du n^gatif utilise 
pour prendre la photographies guelles que soient les 
techniques d fl agrandissement mises en jeu« 
5 Bien entendu les praticiens cherchent Si connaitre 

les images des structures 6tudiees avec la plus grande 
finesse possible s ceci conduit automatiquement a 
realiser, au moment des experimentations, des series de 
mesures avec la plus grande resolution possible «, 

10 L 8 inconvenient presente par une grande resolution et une 
grande finesse reconstruite se situe dans la duree des 
calculs de reconstruction qu°il faut envisages Ceci est 
particulierement le cas avec les methodes itdratives ou 
les calculs sont entrepris plusieurs fois, en general au 

15 moins de deux fois- 

Dans 1 9 invention on resout ce problfeme de la duree 
des calculs en remarquant que, pour les premieres 
iterations, il n'est pas necessaire d'effectuer les 
calculs avec une tres grande resolution ni une tres 

20 grande finesse. En effet, les ecarts evalues au cours 
des premieres iterations sont tellement importants qu B il 
n°est pas necessaire de les evaluer trfes pr6cisemento 
Par exemple, on comprend qu'il n'est pas necessaire de 
connaitre jusqu'a la troisieme decimale la valeur d * une 

25 grandeur physique mesuree dans un voxel quand on sait 
que le resultat sur la premiere decimale peut deja lui 
meme etre faux„ L" invention ne porte cependant pas sur 
le caractere flottant de cette precision qui, si on y 
regarde bien ne conduirait pas a une reduction 

3 0 significative des temps de calculo En effet, les 
operations d' iteration comportent classiquement des 
series de multiplications et d 9 additions dont la duree 
elementaire n f est pas liee au nombre de bit a traiter* 
II y a en effet un traitement parallele en permanence « 
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Cette duree est plutot liee au fait que ces operations 
doivent se suivre i une operation suivante utilisant 
comme arguments les resultats d'une operation 
precedente« 

5 Dans 1 T invention , on a remarque qu'il n'etait pas 

necessaire d'exprimer, pour les premieres iterations, le 
volume reconstruit avec la finesse finale mais que, au 
contraire, on pouvait utiliser une finesse moins grande 
que la finesse finale. En meme temps on utilise une 

10 resolution des mesures moins grande que la resolution 
finale o On attribue le qualificatif final & la finesse 
la plus grande dans laquelle le praticien desire ensuite 
voir les images. Pour la resolution des mesures, la 
resolution finale /est celle qui est utilisee en dernier 

15 dans les calculs (elle est egalement celle qui est 
reellement acquise au cours des experimentations) „ 
Autrement dit, au lieu d'essayer d°exprimer des valeurs 
de grandeurs physiques, en chaque voxel de la structure, 
on essaye de representer cette structure par des 

20 macro-voxels, au cours d B iterations prealables« Avant 
d'effectuer les iterations finales, on transforme ces 
macro-voxels, par une operation d D agrandissement (zoom), 
en des voxels distribues selon la finesse finale o Bien 
entendu pour le calcul des grandeurs physiques a 

25 l'endroit des macro-voxels on utilise des mesures selon 
une resolution moins grande „ 

L f invention a done pour objet un procede de 
reconstruction multi-echelle de 1 D image de la structure 
d'un corps dans lequel 

30 - on effectue, avec une premiere resolution donnee, des 
series de mesures d'au moins une propriete physique 
attachee k des Elements de volume, ou voxels, du corps 
avec, pour chaque serie, une configuration particuliere 
d'une caracteristique d'un protocole donne et en 
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changeant, d'une s6rie a une autre , la configuration de 
cette caracteristique de ce protocole, 

- on acquierfc ainsi des series de resultats de mesure, 
chaque serie etant representative du corps et de la 

5 configuration correspondante de cette caractdristique, 

- et on reconstruit it6rativement une image de la 
structure du corps par mise en oeuvre d'un algorithme de 
reconstruction de type iteratif avec lequel, d'une 
iteration a une autre, on deduit plus exactement une 

10 nouvelle image de la structure du corps en effectuant le 
traitement de reconstruction a partir, d B une part, des 
series de mesures, et d 8 autre part, d'une precedente 
image de cette structure, 

- cette image comportant une collection d 8 elements 
15 d 8 image, chaque element d 8 images comportant une premiere 

information representative du lieu d'un voxel dans le 
corps et une deuxi^me information representative de la 
valeur de ladite propriety physique dans ce voxel, ces 
elements d 1 image etant distribues avec une finesse 
20 donn6e dependant de cette resolution donnee, 
caracterise en ce que 

- on utilise pour au mo ins une iteration donnee des 
informations de mesure correspondant a une resolution 
moins grande que ladite resolution donnee, et une image 

25 de la structure selon une finesse moins grande que la 
finesse donnee, et en ce que 

- pour une iteration suivante, on transforme au 
prealable 1 ' image de la structure selon cette finesse 
moins grande en une image de la structure selon une 

3 0 finesse plus grande que cette finesse moins grande, et 
on utilise en correspondance des informations de mesure 
correspondant h une resolution plus grande que cette 
resolution moins grande „ 

D*une maniere preferee, pour passer des resultats 
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de raesure d fl une resolution finale (acquise des le debut) 
a une resolution moins grande / on procede a une moyenneo 
Cette moyenne peut meme etre ponderee par une fonction 
de lissage, par exemple une gaussienne. Par ailleurs, au 
5 moment ou on trans forme 1 8 image de la structure selon la 
finesse moins grande en une image de la structure selon 
une finesse plus grande , de preference on se livre a une 
interpolation trilineaire* On montrera que cette 
technique de moyenne d'une part, et d 1 interpolation 

10 trilineaire d 8 autre part, conduit d°une part a enlever 
des bruits de reconstruction que l'on aurait 
naturellement obtenus avec les autres proc^des de 
reconstruction o Elle conduit d * autre part , et 
paradoxalement, a une representation plus precise que 

15 celle qu'on obtiendrait avec un calcul de reconstruction 
traditionnel dans lesquels les iterations auraient ete 
faites a chaque fois sur la finesse et la resolution la 
plus grande. On attribue ce resultat superieur aux 
operations de moyennage et d 1 interpolation qui ont pour 

20 effet de lisser les images et d'eliminer naturellement 
les bruits o 

En consequence, en plus d fl aller plus vite, la 
methode de l ff invention permet sur certains types 
d'objets (plus lisse) de dormer un meilleur resultat „ 

25 L" invention sera mieux comprise a la lecture de la 

description qui suit et a l'examen des figures qui 
1 1 accompagnent celles-ci ne sont donnees qu'a titre 
indicatif et nullement limitatif de 1 9 invention «, Les 
figures montrent t 

30 - figure la et figure lb s une representation 
schematique d'un exemple d 8 un mode d ff acquisition et d 0 un 
mode de reconstruction d" images pour lequel on utilise 
le procede de reconstruction de 1' invention? 
- figure 2 s la representation schdmatique d f operations 
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de moyennage, dites aussi de sous-echantillonnage, 
conduisant a des resolutions moins grandes des series de 
mesures acquises; 

- figures 3a a 3c s une representation schematique d B un 
5 mode prefere d°agrandissement de la finesse dans 
lesquels les valeurs des voxels contenus dans les 
macro-voxels sont determinees par interpolation 
trilineaires o 

Les figures la et lb montrent dans un exemple de 

10 densitometrie les conditions d B acquisition de series de 
resultats de mesure, et la reconstruction des images qui 
doit correspondre a cette acquisition en volume « 
L' exemple indique ici correspond a une application de 
densitometrie ou 1 8 invention est particuli&rement utile* 

15 Cependant la methode de 1' invention est aussi applicable 
h d B autres techniques s celles indiquees au d6but de 
cette description o La figure la montre un mode 
d 8 acquisition o La figure lb montre un mode de 
reconstruction » Au moment d " une acquisition de 

20 densitometrie, un tube 1 k rayons X emet un rayonnement 
X 2 en direction d ff un corps 3 d Q uh patient dispose sur 
un panneau porte-patient 4. De 1 8 autre cote du corps 3 
par rapport au tube 1 a rayons X est dispose un 
detecteur 5 c Le detecteur 5 est de type deux dimensions „ 

25 ii possede une pluralite de cellules detectrices telle 
que 6 organisees sur une surface, par exemple 
matriciellement, selon deux directions X et Y* Le tube 
1 a rayon X et le detecteur 5, maintenus par ailleurs 
solidaires l'un de 1 5 autre, peuvent tourner autour d'un 

30 axe 7 passant sensiblement par le milieu du patient 3, 
Les positions en rotation de cet ensemble sont reperees 
par un angle a. A chaque position de l g ensemble tube 1 - 
detecteur 5, on peut acquerir dans chaque cellule 6 des 
resultats de mesures radiologiques « A chaque position on 
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acquiert done une serie de resultats de mesure*, 

Le protocole indique ici est done un protocole de 
densitometrie, et la caracteristique qui conditionne les 
series de rdsultats est la valeur de 1° angle a qui 
5 indique 1 0 orientation de 1 9 ensemble tube 1 - detecteur 5 
par rapport au corps 3 <> Des qu°une acquisition a ete 
faite on fait tourner cet ensemble d 11 tin angle 
elementaire, et on recommence 1 "acquisition en obtenant 
ainsi une deuxieme serie de resultats de mesureo Ainsi 

10 de suite on soumet le corps 3 a des irradiations 
orientees selon un ensemble de directions reparties sur 
au mo ins un demi-cercle- On sait, notamment a partir de 
la demande de brevet frangais citee plus haut, qu°une 
telle serie d ff acquisitions permet de reconstruire une 

I 5 image representative de ce qui se passe dans le corps 3 
a l'endroit ou les rayons X 2 l D ont traverse <, 

La figure lb represente 1 11 algorithme de 
reconstruction par iteration qui est mis en oeuvre dans 
1 5 invention o Dans un tel algorithme on forme l'hypothese 

20 d D une connaissance a priori de 1 D image de la structure 
du corps 3. Cette image 8 comporte un ensemble de voxels 
charges chacun de la valeur d'une grandeur physique o 
Dans cet algorithme , 1 D image 8 est soumise g par 
simulation, a une irradiation provenant d B un point focal 

25 9 situe a une distance compareible a la distance qui 
separe la position du foyer du tube 1 a rayon X par 
rapport au corps 3« Soumise a cette irradiation 
theorique 1° image 8 a trois dimensions de la structure 
se projette en une image 10 a deux dimensions sur un 

3 0 plan de detection occupant une position en 
correspondance avec celle du detecteur 5 par rapport au 
corps 3. Le detecteur simule lO comporte des elements de 
detection comparables aux cellules 6 du detecteur 5 a On 
compare, serie par serie, et, dans chaque serie, cellule 
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a Element, les valeurs des grandeurs physiques 
impliquSes d'une part par les v;iaesures reellement 
acgulses et d* autre part les repro j ections simuleeso On 
en deduit une nouvelle image simulSe 11 a trois 
5 dimensions o Celle-ci est alors substitute h l'ancienne 
image 8 la structure du corps 3o Autrement dit 1' image 
11 de la structure du corps 3 est elaboree a partir de 
la connaissance de 1 ■ ancienne image 8 de la structure de 
ce corps 3, et de la valeur de l'ecart entre la 

10 projection theorique 10 de cette structure 8 et les 
resultats de detection dans les cellules 6, Ceci est le 
principe des algorithmes de reconstruction algebriques 
tels que celui de la demande de brevet citee ci-dessous,, 
On con?oit facilement que la finesse avec laquelle 

15 on est susceptible de decrire 1" image 8 de la structure 
du corps 3 est totalement dSpendante de la resolution 
des elements de detection de la surface lO, et aussi de 
la resolution des cellules sur le detecteur 5. Dans un 
exemple ou le detecteur 5 est un detecteur comportant 

20 256 x 256 cellules on peut montrer qu'on peut 
reconstruire, avec une finesse significative, la 
structure 3 selon un volume avec 256 x 256 x 256 voxels , 
si on a par ailleurs effectu6 256 acquisitions <> 
Normalement, au moment des iterations , les comparaisons 

25 et les calculs de reconstruction doivent etre faits 
pour chacun de ces 256 3 voxels „ 

Dans 1 8 invention , plutot que d'effectuer un calcul 
sur un aussi grand riombre de voxels on procede a une 
simplification., Le principe de cette simplification est 

3 0 represents sur la figure 2. Sur la figure 2 on a 
schematiquement represents, en haut a gauche, la 
repartition des cellules 6 de detection du detecteur 5« 
Ces cellules correspondent aux elements d 1 images de 
detection du detecteur simule 10 o Ces cellules et ces 



2656129 



11 

elements sont distribues avec vine resolution, dans un 
exemple, 256 x 256- 

Dans 1 B invention on decide de reconstruire une 
image de la structure 3 dont la finesse, au lieu d'etre 
5 de 256 3 , ne sera que de 128 3 o Dans ce cas on utilise en 
correspondance une resolution- moins grande des mesures<> 
Dans cet exemple la resolution moins grande est meme de 
moitieo Dans ce but on definit alors des macro-cellules 
de detection, comme la macro-cellule 12 (en haut a 

10 droite) o A cette macro-cellule 12 du detecteur 5 moins 
bien r6solu, on attribue, pour chacune des 256 
orientations de l ff ensemble tube 1 - d6tecteur 5 par 
rapport au corps 3, un macro-resultat de mesure, Ce 
macro-resultat est une combinaison des r6sultats de 

15 mesure dans les cellules 13 a 16 contenues dans la 
macro-cellule 12, Le fait de faire cette moyenne conduit 
ainsi a un lissage dont on a deja indique les effets 
benefiques ci-avanto II est neanmoins possible de 
def inir un voisinage plus grand pour la determination du 

20 macro-resultat attribue a une macro-cellule „ On peut 
prendre par exemple les resultats de mesure dans un 
ensemble de 9 cellules du detecteur 5, comme c°est le 
cas pour la macro-cellule 17 o Dans le cas ou le nombre 
de cellules du voisinage sur le detecteur 5 est impair, 

25 le centre de la macro-cellule sera en correspondance 
exacte avec celui de la cellule centrale 18 autour 
duquel le voisinage a ete etablio Ceci peut conduire par 
ailleurs a des simplifications de balculo II est par 
ailleurs possible de choisir un voisinage de 

3 0 determination du macro-resultat (sur seize cellules par 
exemple) different de l l ensemble (13 a 16) des cellules 
definissant une macro-cellule „ 

Les macro-cellules 12 et 17, quel que soit leur 
mode de combinaison, sont de toutes fa9ons, h un niveau 



2656129 



12 

de calcul donne toutes calculees de la meme fa?on 0 

On a indique precedemment qu'il pouvait etre utile 
d'effectuer au mo ins deux iterations*, Compte tenu de ce 
qu°on cherche une connaissance finale de I 1 image avec 
5 une finesse de 256 3 , on peut choisir de realiser une 
situation avec une resolution encore moins grande, avec 
64 macro-macro-elements ou encore 64 macro-macro- 
cellules au lieu de 128 macro-elements ou 
macro-cellules- Au niveau 64, un macro-macro-element 

10 possede corame information la moyenne des resultats de 
mesures des grandeurs physiques dans un voisinage encore 
plus grand o Par exemple dans le macro-macro-element 19 
(en bas a gauche) du niveau de resolution 64 on a une 
combinaison des resultats dans 16 cellules 20 a 35, 

15 Dans un exemple simple la combinaison des 16 

resultats des cellules 20 a 35 est une moyenne • Dans ce 
cas on peut meme dire que la valeur attribuee au 
macro-macro-element 19 de niveau 64 est la moyenne de 
quatre macro -elements du niveau 128 s les elements 3 6 a 

20 39 0 Dans ce cas 1 * operation de sous 6chantillonnage & 
laquelle il faut soumettre les mesures, pour diminuer 
deux fois la finesse et la resolution, peuvent comporter 
les moyennes sur quatre elements -cellules du d6tecteur 
5, puis les moyennes sur quatre macro-elements-cellules 

25 voisinso 

Cependant, on peut effectuer, pour le calcul de la 
valeur a attribuer aux macro-macro-elements 19 du niveau 
64, une ponderation, par une fonction de lissage, par 
exemple gaussienne, des resultats de mesure des 
3 0 grandeurs physiques dans les cellules 20 a 35. Cette 
ponderation peut alors etre du genre s 

(l/727ra2) 0 exp (- r 2 /2a 2 ) 
Dans cette expression r repr^sente la distance 
entre le centre du macro-macro-element 19 et chacune des 
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cellules 20 a 35 son voisinageo La valeur a est une 
valeur d fl 6cart type fixee comme parametre* Une bonne 
valeur de a est telle que la largeur k mi-hauteur de la 
gaussienne est egale a la taille de la cellule o 
5 A ce stade de 1 8 explication on peut admettre qu ? on 

connait trois groupes de series de mesure„ Dans un 
premier groupe, on a toutes les series de mesures 
finales (celles qui sont ef f ectivement mesurees au 
debut) o Le . nombre de ces series est egal au nombre 

10 d 1 angles a exploration de corps 3« Dans un deuxieme 
groupe, on a compose les mesures de maniere a ce 
qu fl elles se repartissent selon les macro-61ements au 
niveau 128 „ On peut avoir la aussi un meme nombre de 
series o Dans le troisi&me groupe on a, en correspondance 

15 avec les deux premiers groupes, des memes series de 
mesure mais pour des macro-macro-elements du niveau 64 «, 
D'un groupe a 1 5 autre la resolution est de xaoins en 
moins grande (de plus en plus grossiere) o 

Dans 1 fl invention on commence la reconstruction en 

20 effectuant une iteration avec les 64 2 x 256 
macro-macro-elements 19 du niveau 64 o Dans ce cas on 
obtient une image avec une finesse moins grande 
comportant 64 3 macro-macro-voxels „ On utilise ainsi 
toutes les 256 images en projection acquises« II est 

25 neanmoins possible d'en prendre un plus petit nombre o 
Par exemple, on peut en prendre une sur deux (au niveau 
128) ou meme une sur quatre (au niveau 64) » 

On aurait pu calculer 1 8 image au niveau 64 sans se 
livrer aux operations de combinaison et de moyennage 

30 indiquees ci-dessus. Cependant si on agit ainsi, on va 
perdre les petits details, et on risque de voir 
apparaitre dans 1* image finale un tram6 d° image qui 
apparait comme un bruit . 

Une fois qu'une iteration, ou mdme plusieurs 



2656129 



14 

iterations , ont ete entreprises au niveau 64 on peut 
decider d°effectuer une iteration suivante au niveau 
128 o Dans ce cas il est necessaire de transformer 
1° image de la structure 3 avec une finesse en 64 3 
5 macro-macro-voxels en une image avec une finesse en 12 8^ 
macro-voxels o Une solution simple consisterait a 
attribuer a chaque macro-voxel dans une finesse 
superieure les valeurs retenues, pour les 
macro-macro-voxels du niveau inferieur. Cependant f comme 

10 pour les operations de moyennage cette pratique conduit 
a des artefact dans la representation «, Aussi, de 
preference on se livre sur les macro-macro-voxels a une 
interpolation, par exemple trilin^aire, au voisinage de 
cliacun des macro-voxels pour en determiner la valeur., 

15 Les figures 3a a 3c permettent de comprendre le 

calcul de 1 B interpolation trilineaire« La figure 3a 
montre, pour simplifier, le principe du calcul d'une 
interpolation bilineaire. Dans celui-ci, a partir de la 
connaissance de valeurs en des elements d° images A, B, 

2 0 c, D, correspondant a une representation grossiere on 
cherche a connaitre des valeurs a attribuer a des 
sous-elements d fl image, marques par des croix, et plus 
finement distribues„ Pour chacun de ces sous-el6ments 
d 9 images, par exemple le sous-element 40, on mesure ses 

25 coordonnees u, v, dans un repere centre sur 1 Q un des 
elements d" image, par exemple 1° Element A» On calcule 
alors le terme s 

a- (1-u) . (1-v) 

Dans cette expression a est la valeur de 1° element 
30 d" image A, On calcule une meme expression concernant le 
sous-element 40, mais rapportee & chacun des autres 
Elements B, C et Do On additionne ces quatre termes pour 
obtenir la valeur d ■ image du sous-element 40, On 
remarque que ces calculs se simplifient puisque u et v 
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ne peuvent prendre que des valeurs egales a 1/4 et 3/4 
de la distance AC separant deux elements d 5 image entre 
euxo 

La figure 3b montre la connaissance que l'on a de 
5 l 0 image de la structure du corps 3 apres une iteration a 
un niveau inferieur, par exemple au niveau 64« Elle 
montre les centres A a H de huit macro-macro-elements 
d 8 image ou voxel de cette image. La figure 3 c montre la 
connaissance que l g on voudrait avoir, en des 

10 macro-voxels intermediaires tels que 41, de 1 8 image de 
la structure du corps 3 a un niveau de finesse 
super ieure, par exemple 128 a celui du niveau 
precedent o L » interpolation trilineaire dont il est 
question utilise les coordonnees u, v, w des centres de 

15 chacun des macro-61ements d u images 41 par rapport a des 
systfemes d 8 axes centres sur chacun des 

macro-macro-elements d' image au niveau de la finesse 
infer ieure, par exemple de A. La valeur a attribuer au 
macro-element d 5 image 41 est egale alors a la some des 

20 termes tels que % 

a* (1-u) o (1-v) o (1-w) 
dans laquelle a represente la valeur determin6e au cours 
des iterations precedentes pour le macro-macro-voxel A, 
On remarque que dans l 1 invention on a remplace des 
25 multiples calculs d 9 iteration, avec une grande finesse 
et avec une grande resolution, par le calcul 
d B iterations a 1 8 issue desquelles la finesse et la 
resolution augmentent, et aussi par des calculs de 
moyennes d'une part et des calculs d fl interpolations 
3 0 trilin6aires d 8 autre parte On va maintenant montrer que 
le temps de calcul de ces moyennes et de ces 
interpolations est n6gligeable„ 

On sait par exemple que, pour un voxel avec les 
algorithmes de type algebrique, et a chaque iteration, 



V 



2S5S129 



16 

il y a environ 40 operations £t effectuer par projection . 
Le nombre total des operations est done, pour une 
iteration donnee, directement proportionnel au nombre de 
voxels 64 3 , 128 3 , ou 256 3 . Dans un exemple, pour fixer 
5 les id6es, on remplace une iteration a 128 suivie d'une 
iteration a 128 , par une iteration a 64 suivie d'un 
agrandissement de 64 a 128 et d'une iteration a 128. Si 
on appelle n p le nombre de projections dans la serie, et 
n v le nombre de voxels f on peut dire qu'une iteration a 

10 un niveau donne (par exemple 64) conduit a un nombre 
d 5 operations egal a 2 x 40 x n p xn v « Par contre, 
1 8 iteration a un niveau superieur (par exemple en 128) 
conduit a un nombre d 11 operations egal a 2 x 40 X np x n v 
x 8o En effet etant donne qu'il y a huit fois plus de 

15 voxels en 128 qu D il n"y en a en 64 f il y a huit fois 
plus de calculs<, 

En simplifiant on peut admettre qu^l y a en 
moyenne moins de 6 operations d l addition ou de 
multiplication pour determiner par interpolation 

20 trilineaire, la valeur en tin macro-element d' images 41 
quand on connait la valeur en 8 macro -macro-elements 
d' image voisins A a Ho Le nombre total d 8 operations 
necessaires pour faire 1 ff agrandissement est alors egal a 
6 x n v x 8. En effet, ce nombre est proportionnel k huit 

25 fois le nombre de macro-macro-voxels „ Un rapide calcul 
montre qu'une iteration a 64 suivie d'un agrandissement 
donne un nombre d 1 operation tout a fait negligeable par- 
rapport au nombre d 8 operations avec deux iterations en 
128 o 

30 Autrement dit dans 1° invention on fait f pendant 

une duree de calcul sensiblement egale a une iteration a 
128, une iteration a 64 et une iteration a 128 « On 
montre bien ainsi qu'un rapport sensiblement 2 est gagne 
sur le temps. En pratique il y a une premiere iteration 



2656129 



17 

en 64 pour determiner une image en finesse 64 3 o Par 
interpolation on deduit une image en finesse 128 3 o On 
effectuera sur cette image 128 3 une iteration h 128 0 
Apres l 1 iteration en 128 , on deduit une image en finesse 
5 256 3 par agrandissement. Sur cette image en 256 3 on 
effectue la derniere iteration „ Le temps de calcul est 
alors celui de la derniere iteration seule« Le temps 
total est ainsi fortement reduit pour une meme qualite 
d ■ image * 

10 si on ne procedait pas par interpolation 

trilineaire on pourrait aussi affecter & chaque 
macro-voxel 41 , appartenant au voisinage d ff un 
macro-macro-voxel A r la valeur a retenue pour A„ Dans ce 
cas cependant on observe un phenomene de pavage entre 

15 les macro-voxels du voisinage de A et ceux du voisinage 
de B„ Ce phenomene apparait comme V tin artefact dans 
1 8 image d ff une maniere similaire au sous echantillonnage 
sans la moyenneo Une autre solution que 1 1 interpolation 
trilineaire (par exemple ponderee par une gaussienne) 

20 est la aussi envisageableo 

Les m6thodes de reconstruction algebriques 
pr6sentent par ailleurs un autre avantage« Elles 
permettent de determiner un support de l fl objet a 
reconstruire o Cette determination a essentiellement pour 

25 objet de limiter le calcul de 1 ? image de la structure de 
l B objet en des elements d 1 images qui correspondent 
seulement a des endroits ou se trouvent des parties 
interessantes de cet objet „ Par exemple, dans le corps 
humain, et en angiographies on ne reconstruit que 

3 0 l 8 arbre vasculaire, pas les tissus qui l B entourent« Dans 
ce but on commence par definir un support de l'objeto 
En effet, ces structures sont des structures "creuses" 
et il est inutile de calculer les valeurs des structures 
environnantes , pour les eliminer ensuite au cours d g une 
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operation de visualisation de la structure reconstruite „ 
On sait, notamment par la demande du brevet frangais 
n°88 17001 depose le 22 decembre 1988 par le meme 
demandeur, definir rapidement une region support d 8 un 
5 objeto 

La question se pose alors dans 1 0 invention de 
savoir comment va etre defini le support , compte tenu du 
fait que dans les premieres iterations on utilise une 
finesse moins grande de 1 8 image 8 de la structure 3„ On 

10 s'est alors rendu compte qu'il etait necessaire, pour 
definir le support pour une petite finesse, d°utiliser 
la connaissance du support pour une grande finesse * En 
consequence on doit considerer comme appartenant au 
support, a une finesse moins grande (128 3 ) , tous les 

15 macro-voxels qui contiendraient ne serait-ce qu g un voxel 
de 1" image definie sur la grande finesse (256 3 ) <> En 
effet, si on calculait specif iquement un support au 
niveau de chaque finesse (64 ou 128 ) on risquerait 
d°eliminer, pour les premieres iterations, des 

2 0 contributions d 8 elements d 8 image que l B on ne pourrait 

jamais incorporer correctement dans le calcul par la 
suite « 

Par ailleurs, on s 8 est rendu compte qu°il dtait 
tres int6ressant d'effectuer des traitements d° image sur 
25 les images intermediaires (selon une finesse moins 
grande) <, En effet, lorsqu'on veut utiliser une image on 
cherche, en particulier pour la visual iser, a en enlever 
le bruit o Les differentes techniques de traitement qui 
permettent cette amelioration consistent a filtrer les 

3 0 elements d 8 image ou encore a segmenter I 9 image « Avec le 

filtrage on attribue comme valeur de grandeur physique a 
un element de volume une valeur dependant de sa propre 
valeur et aussi des valeurs propres attributes aux 
voisins de cet element de volume <, Avec la segmentation 
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on elimine de la structure les elements de volume dont 
la valeur de la grandeur physique est trop faible 
(s°apparente a du bruit) ou, autrement, appartient a une 
plage de valeur inint6ressante (elimination des tissus 
5 mous dans certaines representations osseuses par 
example) 0 

On distingue, parmi les filtrages, d'une part les 
filtrages lindaires, d' autre part, les filtrages non 
lineaires* Par exemple, par un filtrage lir\eaire par 

10 lissage ou moyennage on peut favoriser la reconstruction 
des structures lisses comportant des parties avec une 
densite qui varie lentement ou meme qui est constanteo 
Dans un autre exemple, avec un filtrage non lineaire, on 
va au contraire favoriser la reconstruction des objets 

15 dans la structure dont les valeurs de grandeur physique 
n ff appartiennent pas a des plages de valeur 
interpenetreeso Ces plages de valeur sont sans continuum 
entre elleso 

Quand on segmente, on elimine de la region support 
20 de la structure tous les elements de volume de cette 
structure dont la valeur de la grandeur physique mesuree 
sort d°une plage predeterminee« 

Tous ces traitements, lorsqu°ils sont effectues sur 
les images intermediaires en finesse moins grande, 
25 conduisent k des gains de temps sur le traitement 
lui-meme et sur la reconstruction « Le gain de temps sur 
le traitement est lie au fait qu*on traite moins 
d 8 elements de volume avec une - finesse moins grande 
qu'avec une finesse plus grande « Le gain de temps en 
30 reconstruction est li6 au fait qu ff apr6s segmentation 
intermedia ire il y a moins d°616ments de volume a 
reconstruire (done moins de calculs) que si la 
segmentation est faite sur l'objet reconstruit global- 
On pourrait montrer avec l ff invention que le co<it de 
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temps de ces traitements interm6diaires est largement 
compense par les gains de temps sur les traitements 
finaux. 
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REVENDICATIONS 



5 i - Proc6d6 de reconstruction multi-echelle de 1° image 
de la structure d°un corps dans lequel 

- on effectue, avec une premiere resolution donnee, des 
series de mesures d'au moins une propriete physique 
attachee a des elements de volumes, ou voxels, du corps 

10 avec, pour chague serie, une configuration particuliere 
d B une caracteristique d°un protocole donne et en 
change ant, d'une serie a une autre, la configuration de 
cette caracteristique de ce protocole, 

- on acquiert ainsi des series de resultats de mesure, 
15 chaque sdrie etant representative du corps et de la 

configuration correspondante de cette caracteristique , 

- et on reconstruit iterativement une image de la 
structure du corps par mise en oeuvre d fl un algorithme de 
reconstruction avec lequel, d B une iteration a une autre, 

20 on deduit plus exactement une nouvelle image de la 
structure du corps en effectuant le traitement de 
reconstruction a partir d'une part des series de mesures 
et d° autre part d°une pr6cedente image de cette 
structure, 

25 - cette image comportant une collection d 0 elements 
d° image, chaque element d tt image comportant une premiere 
information representative du lieu d'un voxel dans le 
corps et tine deuxieme information representative de la 
valeur de ladite propriete physique en ce voxel, ces 

30 elements d 8 images 6tant distribues avec une finesse 
donn£e dependant de cette resolution donnee, 
caracterise en ce que 

- on utilise pour au moins une iteration donnee des 
informations de mesure correspondant a une resolution 
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mo ins grande que ladite resolution donn6e, et une image 
de la structure selon une finesse moins grande que la 
finesse donnee, et en ce que, 

- pour une iteration suivante, on trans forme au 
5 prealable I 9 image de la structure selon cette finesse 

moins grande en une image de la structure selon une 
finesse plus grande que cette finesse moins grande, et 
on utilise en correspondance des informations de mesure 
correspondant a une resolution plus grande que cette 
10 resolution moins grande « 

2 - Procede selon la revendication 1, caracterise 
en ce que 

- le protocole est un protocole de type densitom&trie h 
rayons X, 

15 - la caracteristique concerne 1 * or ientation d'un 
densitometre servant a mettre en oeuvre ce protocole 
densitometrique ; 

- la configuration concerne la valeur d°un angle 
d 0 orientation de ce densitometre par rapport a une 

20 reference li6e au corps, 

- les series de mesures comportent des series 
d° acquisitions radiologiques au cours desquelles la 
structure est projetee en une projection & deux 
dimensions sur un d6tecteur a deux dimensions comportant 

25 des cellules d^tectrices reparties selon ladite premifere 
resolution, 

- 1 8 image de la structure du corps est une image & trois 
dimensions, 

- et en ce que la methode de reconstruction est du type 
3 0 algebrique „ 

3 - Procede selon la revendication 2, caracterise 
en ce que, pour la resolution moins grande des 
informations de mesure, on choisit de consid6rer comme 
resultat de mesure, a chaque projection, une combinaison 
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de resultats de mesure en des cellules adjacentes entre 
elles, cette combinaison etant affectee a une 
macro-cellule representative de ces cellules adjacentes <> 

4 - Precede selon la revendication 3, caracterise 
5 en ce que la combinaison est une moyenne arithmetique., 

5 - Procede selon la revendication 4, caracterise 
en ce que la combinaison est une moyenne pondere par une 
fonction de lissage 

6 - Procede selon l'une quelconque des 
10 revendications 1 a 5, caracterise en ce que pour 

transformer 1 8 image de la structure selon une finesse 
mo ins grande en une image selon une finesse plus grande 
on effectue pour des macro-elements (41) selon cette 
finesse plus grande , une interpolation trilindaire de 
15 valeurs (a,b,«<,o) de ladite propriety physique entre des 
macro-macro-elements voisins (A,B, ••-•-H) et distribues 
selon la finesse moins grande 

7 - Procedd selon 1 11 une quelconque des 
revendications 1 a 6, caracteris6 en ce que on limite la 

20 reconstruction iterative a une region support detectee 
en fonction de la finesse plus grande « 

8 - Procede selon l D une quelconque des 
revendications 1 a 7, caracterise en ce que le rapport 
des finesses et des resolutions plus grande et moins 

25 grande sont des multiples de 2<, 

9 - Procede selon l*une quelconque des 
revendications 1 a 8, caracterise en ce que, avant une 
iteration suivant une iteration donnee, on traite 
1 0 image de la structure selon la finesse moins grande de 

3 0 fagon a ameliorer la quality de la reconstruction « 

10 - Procede selon la revendication S, caracteris6 
en ce que pour le traitement, on filtre par lissage les 
elements de volume de la structure selon la finesse 
moins grande de fagon a amdliorer la reconstruction d D un 
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objet physiquement lisse» 

11 - Procede selon la revendication 9, caract^rise 
en ce que pour le traitement, on filtre par lissage non 
lineaire les Elements de volume de la structure selon la 
5 finesse mo ins grande de fa?on a preserver les contours 
et a. ameliorer la reconstruction d ff objets dans la 
structure dont les valeurs des proprietes physiques 
appartiennent a des plages dif ferenci6es et sans partie 
commune - 

10 12 - Procede selon la revendication 9, caracterise 

en ce que pour le traitement, on segmente les Elements 
de volume de la structure selon la finesse mo ins grande 
de maniere k eliminer d 8 une region support de cette 
structure tous les elements de volumes dont la valeur de 

15 la grandeur physique s 8 apparentent a des bruits « 
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